Exact Ground States of One-Dimensional Quantum Systems: 

Matrix Product Approach 



Gang Su* 

Institut fur Theoretische Physik, Universitat zu Koln 
Ziilpicher Strasse 77, D-50937 Koln, Germany 



Abstract 

By using the so-called matrix-product ground state approach, a few one- 
dimensional quantum systems, including a frustrated spin-1/2 Heisenberg 
ladder, the ferromagnetic t-J-V model at half- filling, the antiferromagnetic 
J z — V at 2/3 filling and the antiferromagnetic t — J z — V model at half- 
filling, are solved exactly. The correlation functions in the ground states are 

calculated respectively. Some relevant results are also discussed. 
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I. INTRODUCTION 



Recently there has been growing interest on one- dimensional (ID) quantum systems due 
to a variety of reasons. Since the difficulty in dealing with many-body problems the methods 
in obtaining exact results in ID are rare. As is well-known, a few ID many-body problems 
can be exactly solved by the Bethe ansatz, as the systems satisfy the so-called Yang-Baxter 
equation ffl. Actually, Bethe ansatz is a powerful and efficient method in obtaining exact 
results including ground state and excited states as well as thermodynamics. However, a lot 
of ID quantum systems does not obey Yang-Baxter equation, and thus are non-integrable. 
In this situation, it is difficult, due to absence of a systematic method, to obtain some 
exact results of ID many-body systems. Amongst very few existing exact methods there 
is a so-called matrix product (MP) ground state approach which can be applied to obtain 
exact ground-state properties of ID quantum systems. This method was established a few 
years ago by Kliimper et al where they studied the exact ground state of a large class of 
antiferromagnetic spin-1 chains. In those cases the method was really efficient. The basic 
idea of the method is that the global ground state can be constructed in the form of a 
matrix product of single-site states. This state is usually said optimal in the sense that it 
is the product of all the local ground states of the local interaction hjj+i (see below). To 
show that the MP state is really the ground state, one can diagonalize hjj+x exactly in all 
local eigenstates. If the eigenvalue of hjj+i in the MP ground state is not larger than the 
minimum of all eigenvalues of hjj+i diagonalized in all local eigenstates, then the MP state 
is the global ground state of the whole Hamiltonian. The similar spirit has been applied to 
construct the exact ground states of generalized Hubbard models 0. In this paper, we will 
apply the MP approach to a few quantum systems including a frustrated spin-1 /2 Heisenberg 
spin ladder, and several simplified cases of t-J-V model. The correlation functions of the 
ground states are obtained. Some relevant results are also addressed. 
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II. A FRUSTRATED SPIN-1/2 HEISENBERG LADDER 



Investigations on spin ladders have attracted a lot of attention in recent years 0- 12 



This is usually motivated from the following possible reasons: one is that the Heisenberg 
spin ladder (two antiferromagnetically coupled Heisenberg antiferromagnetic spin- 1/2 chains 



or two ferromagnetically coupled such chains provides a simple model to gain insight 

into the underlying physics of crossover from ID chain to 2D square lattice; one is due to 
the Haldane's conjecture [I5| — which states that a quantum antiferromagnet with integral 
spin has different properties from one with half- integral spin, i.e., the former has a unique 
ground state, a gap to the excited state and exponential decaying ground-state correlation 
functions, which have been confirmed in experiments, while the latter does not — people 
believe that the Heisenberg spin ladder is probably in the same phase as the Heisenberg 
spin-1 antiferromagnetic chain |f),l0|,12|; and the other is that spin ladders are intimately 
related to magnetic properties of some realistic materials (for example, Sr n -iCu n+ i0 2n 
|16|,|8[| and (yO^iPiO-j fTTH). Besides, frustrated spin ladders are also interesting, because 
they are believed to be responsible for physical properties of stoichiometric Sr n ~\Cu n +\Oi n 
compounds ||. Exact solutions on spin ladders are rare so far, whereas earlier studies on 
these models are either approximate (mean-field treatment or bosonization) or numerical, 
although some consensus is qualitatively made. Any exact result on these models will 
hence be necessary. In this section, we report that a frustrated spin-1 /2 Heisenberg ladder, 
with properly constricted parameters, can be solved exactly by using the MP ground state 
approach. 

We consider a spin-1/2 Heisenberg ladder with diagonal coupling, described by the fol- 
lowing Hamiltonian 

H = J i E E s»0') • s nti + 1) + -h E SxOO • s 2 (j) 

71=1,2 j=l j = l 

+j s E[Si(j) ■ s 2(j + 1) + s 2 (j) ■ SiO" + 1)] 

3=1 
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| 2(J 1 + J 3 ) + J^ (1) 

where j runs over all rungs and n over two legs: 1 and 2, L is the length of the ladder, S n (j) is 
spin-1/2 operator on the rung j along leg n of the ladder, J\ is the interaction along the legs, 
J 2 is the coupling along the rungs, and J 3 is the diagonal coupling between the two nearest 
neighbor rungs. A constant (the last term of (1)) is added for convenience of our purpose. 
The periodic boundary condition is assumed along the ladder. To guarantee the solvability 
of the ladder, we have to confine the parameters to satisfy the following conditions: 

0<Ji<- T J f' , J 2 <0, J 3 <0. (2) 

Under the conditions of (2), one may see that spins on the ladder are frustratedly distributed 
along the legs and rungs, minimizing the total energy of the system. We mean the frustration 
in the present model just in this sense. In the following we for brevity denote the local 
eigenstates | + \)j and | — \)j of S z (j) simply by and \—)j, respectively. We denote a 
local state on each rung by \a)' l J l \b) r - , where {a, b} = {+, — }, and {m, n} = {1,2} identifying 
the legs. Introduce a permutation operator 



P^ n > = -[1 + AS m (j) ■ S n (l)]. (3) 
The operator has the property: 

Pfr n) \a)TW = W\a)?. (4) 
With the aid of (3), we rewrite (1) as 



L 



tf = E*«+i> ( 5 ) 

3=1 



with the local interaction 



h - Jl (p( 1 > 1 ) _|_ p( 2 > 2 h , J 2,p(l>2) ( i j2) J 3 (p (l,2) p (2,l)x 

where we have made use of the periodic boundary condition along the ladder. 



Following the discussions in Ref. 0, and using the local eigenstates and \ of S z 



j ■ 



we define a local 2x2 matrix on each rung by 

(7) 



^ !+>}!->? + !->}!+)? l+>}|+>? + l-»l->? ^ 



v !+>}!+>? + l->}l->? !+>}!->? + !->}!+>? y 



9j 

and a global state 

|*o) = Tr gi ®g 2 ® •••®# L (8) 

where ® denotes usual matrix multiplication of 2 x 2 matrices with a tensor product of the 
matrix elements. It can be found that 

ff|*o) = £o|*o>, (9) 
E = ( Jx + y + J z) L - ( 10 ) 

As we have the constraint (2), thus E < 0. Now let us show |\P ) is a ground state of H. 
To realize this purpose, we in turn need to prove gj <S> Qj+i and thus |^ ) to be the ground 
state of hjj+i. After a little tedious algebraic calculations, one may find that 

hj, j+ i(gj ® g j+ i) = ( Ji + y + J 3 )(& ® 

namely, ^ ® Sj+i is an eigenstate of fojj+i with eigenvalue E /L. On the other hand, hjj+i 
can be exactly diagonalized in all possible 16 eigenstates on a plaquette consisting of two 
rungs (j and j + 1) and two corresponding legs. It is not hard to obtain the eigenvalues 
as follows: J x (3-fold), f (3-fold), J 3 (3-fold), J x + f + J 3 (5-fold), ±[J* - § J x J 2 + ± J| - 
-^1-^3 — |<^2^3 + Jf] 5 - One can easily check that E /L = J\ + 4jp + J 3 is the lowest eigenvalue 
between them, under the conditions of (2). Thus, ^ ® is really the ground state of 
hjj + \, and therefore |\I/o) is also the ground state of hj 7 j + \, which in turn proves that |^o) 
is the global ground state of H, and E is the ground-state energy. 

Similar to the calculations in Ref. ||, we can obtain correlation functions in |\I/o) state 
by using the transfer matrix method. The norm (\J/q I ^o) can be obtained by 
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(^dl^o) ^ 9 ni n29n2Ti 3 ' ' ' 9n L n 1 9m 1 m29m2m 3 ' ' ' 9m L mi TrG 



(11) 



{n a ,m a } 

where G is the transfer matrix, with elements G Ml/Lt2 = G ? ( ni . mi ) i ( ri2rrt2 ) = fl , ^ in2 5 , mim 2 where 
the ordering of multi-indices is chosen as usual: \i — 1, 2, 3, 4 *-> (11), (12), (21), (22). The 
4x4 transfer matrix has four eigenvalues: (2- fold) and 4 (2- fold). Thus, one can write 
G\e n ) = A n |e n ), with A n (n = 1,2,3,4) the eigenvalues (Ai j2 = 0, A 3)4 = 4), and \e n ) the 
corresponding normalized eigenvectors which have alternative forms 



/ i \ 

2 



\ei) 



( _A 

2 



|e 2 ) 



2 



|e 3 ) 



V 2 J 



2 



l e 4) 



(12) 



One may verify that {|e n )} forms a complete set, which can be used to calculate the trace 
in the following. Evidently, ($o|^o) = 2 • 4 L . The one-site expectations can be evaluated by 



(A) = ( ? T ol ff v 0> = {TrG L )- l TrZ{A)G L -\ 



(13) 



(*o|*o) 

with Z(i) W)J2 = <?™ n2 ^4<?mim 2 - The 2-site correlations of operators A(l) at site 1 and B(r) 
at site r can be calculated by 



(A(l)B(r)) = (TrG L )- 1 TrZ(A)G r - 2 Z(B)G 



iL-r 



(14) 



In accordance with (12)-(14), we can obtain the following expectations: 



(s^U)) = o, (s z tot ) = E(^O') + sm = o, (s+ 2 (j)) = o, 

3=1 

(SU 1 )^)) = (S z 1>2 (l)SI A (r)) = 0, 

((SUm = (St, 2 (m, 2 (r)) = (S+ 2 (l)S^(r)) = ], 



(15) 
(16) 
(17) 



for r > 2, where the technical details refer to Ref. [[|. We observe that in the ground state the 
model has transverse fluctuations for the transverse correlation functions are uniform with 
varying spatial distances, implying the correlation length for transverse magnetic orderings 
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is infinite, while it does not exhibit longitudinal magnetic order and fluctuations. As can be 
seen, |\& ) is XY magnetic ordered. 



To identify the structure of the state, one can use a matrix u = 

V -v 



to make 



a similarity transformation on gj as g', = ugju 1 



( a + b 
y a — b j 



with a = l+yA-) 2 , + 



|-))|+)| and b = |+))|+)| + Thus |^ ) = Trg[ ® g' 2 <g> ■ ■ ■ ® g' L = \a + 6>i|a + 

b) 2 ■ ■ ■ \a + b) L + \a — b)i\a — b)2 ■ ■ ■ \a — b) l- Obviously, the global state is a linear combination 
of a direct product of all possible spin configurations on each rung, and thus has transverse 
fluctuations but without longitudinal ones, as revealed by the expectations. In addition, 
since gj is a mixture of singlet and triplet states, the excitation from the ground state |^ ) 
may be gapless. 

A simple argument can show that the ground state is degenerate. For instance, if one 
replaces the four entries of matrix gj in (7) by the same expression: or |— )|, one 

will find that the replaced MPG state has the same eigenvalue as that of |\& )- Actually, such 
states are nothing but the fully polarized ferromagnetic states, which can be implemented 
through a similarity matrix transformation on gj. One can also compute expectations of 
some physical quantities in the polarized ferromagnetic state, and will find that they have 
different values as in |^ ) state, while the longitudinal correlation functions are still uniform 
in spatial space. Hence |\& ) an d the fully polarized ferromagnetic state are not the same, 
suggesting that the ground state is not unique. 

That the model has the fully polarized ferromagnetic ground state is not surprising, 
because the coupling along the rungs ( J 2 < 0) and the diagonal coupling between the rungs 
(J3 < 0) favour to form triplets and thus prefer to form ferromagnetic order, while the 
strength along the legs, though it is positive and favours to form antiferromagnetic order, 
is smaller than the other two interactions. The competing result would lead to favour 
ferromagnetism. The present result just confirms the intuitive argument. However, the 
XY magnetic ordered ground state exist in the model is yet plausible in physics. Though 
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the spins on each rung of the ladder prefer to form triplets, minimizing the energy of the 
system, the alignment of the spins along positive direction has the equal weight as that along 
negative direction, and thus the resultant state may have transverse magnetic fluctuations 
but without longitudinal magnetic order, as can be seen clearly from the explicit form of |^o)- 
As the ground state is degenerate, one may anticipate that the spin rotational symmetry is 
spontaneously broken in the ground state, as it should be. To this end, one may note that 
the advantage of MP ground state approach applying to this model lies in having found a 
ground state with transverse magnetic fluctuations, thereby showing the ground state of the 
model is not unique, and the excitation is gapless. 

To conclude this section, the frustrated spin- 1/2 Heisenberg ladder, defined in (1), is 
solved exactly by using the MP ground state approach. It is shown that the ground state of 
the model is degenerate, as a magnetic ordered state with transverse fluctuations and a fully 
polarized ferromagnetic state coexist in the ground state. The excitation from the ground 
state may be gapless. The correlation functions are uniform in spatial space, implying that 
the correlation length is infinite, and thus the system is critical in the ground state. One 
may observe that this model thus provides an example that quasi one-dimensional frustrated 
Heisenberg model can exhibit magnetic order in the ground state, although the Mermin- 
Wagner thermodynamic fluctuations enable the system probably not to possess such an 



order at finite temperatures. We expect that the present result, on the one hand, may shed 
some light on the theoretical basis of some quasi one-dimensional ferromagnets, and on the 
other hand, may provide a clue to experimentalists to synthesize new quasi one- dimensional 
ferromagnetic materials, particularly organic magnets. 

III. T-J-V MODEL 



This model, as an extension of the usual t-J model, first proposed by Schlottman ||T9[ for 



heavy- fermions, and later studied numerically by others [E(T-|22[ . Apart from the supersym- 



metric points |T9[, this model is non-integrable. The system is described by the following 



S 



Hamiltonian 

L 



H tJV = J2 h 3,j+ii ( 18 ) 
i=i 



with the local interaction 



a & 

on ID chain with the length L even, where Cj t<T is the annihilation operator for an electron 
with spin a (=|, I) at site j, {Sj} are spin-1/2 operators, rij = J2a c j,o- c j,a is the electron 
number operator at site j, t is the hopping matrix element, J, J z are anisotropic spin- 
exchange interactions, and V is the strength of density-density interaction. The periodic 
boundary conditions are assumed. We further suppose that, as usual, the double-occupancy 
of every site is forbidden due to the existence of a large on-site Coulomb repulsion, i.e., 
each site has either one electron (with spin up or down) or empty. Therefore, there are 
three possible electronic states at a given site j: |0), the Fock vacuum satisfying c Jj0 - J 0) = 0; 
| |) = cj JO); and | |) = cjj ,|0). Since there exists such a property, the problem, as first 



noted by Schlottman [F9| in the t-J-V model, can be rewritten in terms of spin operators 
corresponding to S = 1 in the restricted Hilbert space, which implies that the model should 
have similar behaviors as that of S = 1 spin chain. This defines the t-J-V model. Since 
it is non-integrable off supersymmetric points, it is difficult to obtain the exact results of 
this model. In the following we will concentrate on a few special cases, which can be solved 
exactly using MP approach. 



A. Ferromagnetic t-J-V model at half-filling 



In this case we assume that J z = J < and t > 0. Similar to Sec. II, we define a 2 x 2 
matrix 



' IT) + 11) ^(1 T>- 

v 7^(11) -IT)) ID + IT) j 



(20) 



with a > a parameter, and propose an ansatz for the global ground state 

\$i)=Trg , 1 ®g , 2 ®---®g , L , (21) 

where the notation is the same as in Sec. II. Following exactly the analyses in Sec. II, it 
is readily show that is the global ground state, with eigenvalue (V + ^)L, of the 
ferromagnetic t-J-V model under the following constraint: 

t>0, J<0, V<-t-j. (22) 

The expectation values in the ground state can be calculated directly using the transfer 
matrix method. Though the finite-size results are available, we are only interested in those 
in the thermodynamic limit. Henceafter we mean the result just for L — > oo. It can be 
found that 

(S-) = 0, (SZ*) = (£%) = 0, (fif> = I^i, („,.) = !, ((S*f) = \, (J,) = 0, (23) 

with Jj = itJ2a(. cJ j+i,a c j,(T~ c },cr c j+i,o-), the current density. Therefore, we see that the ground 
state has XY ordering and is insulating. The two-point correlation functions are evaluated 

(5?5;>=0, {S+S-) = (Sffi*) = 1 {nwr) = l, (mS;) = 0. (24) 

When a — 1, (Sf) = 0, implying that there is a phase transition from one XY ordered 
state with nonvanishing order parameter into another XY ordered state in which the order 
parameter vanishes while the transverse correlation function is uniform. Since the correlation 
length of transverse correlation functions is infinite, the system is critical in the ground state. 
An interesting observation is that, the state is still the ground state, if J = so long as 
V < —t. This fact suggests that the density-density interaction seems to play a significant 
role in the mechanism of the origin of ferromagnetism. Similar to arguments in Sec. II, the 
ground state of the system may be degenerate with the fully polarized ferromagnetic state. 
Since the symmetry breaking requires the degeneracy of the ground states, the present result 
is reasonable, as (23) shows that the spin rotational symmetry is spontaneously broken in 
the ground state. If one makes a particle-hole transformation in the system, then he will 
find the transformed system will exhibit off-diagonal long-range order in the ground state. 
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B. Antiferromagnetic J z — V model at 2/3 filling 



Define a local 2x2 matrix 

> + m + 10} i n - 1 \) - \o) \ 

(25) 



9" 



Vii)-it) + io) -ii)-it)-|o> 

and an ansatz for the global ground state of H t jv 

\<$> 2 ) = Trg'{®g'l®---®g'l (26) 

To make |<E>2) be the global ground state of H t jv, tracing the exact way in Sec. II, the 
following conditions must be satisfied: 

t = J = 0, V = ^j, J z >0. (27) 

One can see that, |$2) is actually the global ground state, with eigenvalue 0, of antiferromag- 
netic Ising model with density- density interaction described by the following Hamiltonian 

Hj.y = J Z J2 S*S] +1 + ^E n jnj+1 . (28) 

3=1 ^ j=l 

It is well-known that this model without the density-density interaction has been solved by 
Ising [23] seventy years ago. Now we have applied matrix product to construct the global 



ground state of the Ising model with density- density interaction. Applying the similar 
arguments in Ref. one can show that the ground state is unique. One will see later that 
the state is at 2/3 filling. The state is thus nontrivial and is optimal. 
The expectation values in the ground state can be obtained: 

(S!) = 0, (Sl«) = (t,Sj) = 0, {S}) = 0, (n,) = 2 -, ((S]) 2 ) = l. (29) 

The spin-spin correlation function is nontrivial 

(S'S?) = 0.0331456 - 0.0625(-l) r , (r > 2) (30) 

which means that the global ground state is antiferromagnetic, and has long-range order. 
Since the system is localized, the ground state should be insulating. The spin-charge corre- 
lation function is 

11 



(S^n r ) = 0.0331456 - 0.04167(-l) r , (r > 2) (31) 

which suggests that the spin and charge are not separated in the ground state. The other 
correlation functions are: 

(S+S;) = (S^S+) = 0, ( ni n r ) = 0.5, ((^) 2 (^) 2 > = 0.03125. (32) 

These correlation functions are nontrivial, and it seems that they appear for the first time. 
One may note that these properties are quite different from those of Ising model Hlip3|. 

C. Antiferromagnetic t — J z — V model at half- filling 

Consider the simplified version of H t jv, described by the Hamiltonian (18) with J = 0, 
which we call the t — J z — V model. Below we will show that the ground states of this 
simplied model can be constructed using MP ground state approach at half-filling. As 
before, we define an 2 x 2 matrix 



/// 

9n 



'|T) + I0 IT) -I I) ^ 
J I) -IT) -I I) -IT) J. 



(33) 



We propose an ansatz 

|d> 3 > =Trg'{'®g^®---®g'^ (34) 

Pursuing the same procedure as in Sec. II, one can find that |$3) is really the global ground 
state, with eigenvalue (V — ^f)L, of the t — J z — V model, subject to the following constraint 

J z > 0, t > 0, V < -t + (35) 

Similar to the arguments in Ref. 0, it can be shown that the global ground state is unique. 
The expectation values in the ground state are obtained 

W> = o, (S£*> = <!>/> = o, (sf> = o, K-) = i, ((s;) 2 ) = 1 -. (36) 

3=1 4 
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The global ground state is thus antiferromagnetic. The spin-spin correlation function is 
found to be 

<^> = -i(-ir, (r>2) (37) 

which is nothing but the result of Ising model |T]J2^] . This fact shows that the antiferromag- 
netic t — J z — V model at half-filling has the same properties as the Ising model, implying 
that the model at half-filling can be mapped to the Ising model. Actually, this is true, be- 
cause the hopping term, at half-filling, plays no role in the restricted Hilbert space without 
double-occupied sites. The other two-point correlation functions are 

(S+S;) = (S^S+) = 0, (mn r ) = 1, (S*n r ) = 0, ((Sf) 2 ^) 2 ) = 0.03125. (38) 

One may observe that the behaviors of correlators at half-filling are quite different from 
those at 2/3 filling, as discussed in last subsection. 

IV. SUMMARY AND CONCLUSION 

We have studied the exact ground-state properties of a few ID quantum systems including 
a frustrated spin- 1/2 Heisenberg ladder, a ferromagnetic t-J-V model, an antiferromagnetic 
J z — V model at 2/3 filling and an antiferromagnetic t — J z — V model at half-filling, by 
using the matrix product ground state approach. The correlation functions were obtained 
in the ground states. 

It is shown that the ground state of the frustrated spin- 1/2 Heisenberg ladder with 
restricted parameters is magnetic ordered and is degenerate. The excitation from the ground 
state may be gapless. The correlation functions are uniform in spatial space, implying that 
the system is critical in the ground state in the sense that the correlation length is infinite. 

The exact ground state of the ferromagnetic t-J-V model with proper parameters was 
constructed at half-filling. It is found that the ground state has magnetic long-range order, 
and is also degenerate. A phase transition is found in the system. The spin rotational 
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symmetry is spontaneously broken in the state. The correlation function is uniform, implies 
that the system is also critical in the ground state. 

The global ground state of the antiferromagnetic J z — V model with V = J z /4: is found 
to be antiferromagnetic and unique at 2/3 filling. The spin-spin correlation functions is fluc- 
tuating, non-vanishing eventually, with increasing the spatial distances, suggesting that the 
ground state has antiferromagnetic long-range order. The spin-charge correlation function is 
also fluctuating, implying that the spin and charge degrees of freedom are highly correlated 
in the state. The spin and charge are not seperated. The behaviors are quite different from 
those of Ising model. 

The studies on the antiferromagnetic t — J z — V model with proper parameters at half- 
filling show that the system is very similar to the Ising model, as the spin-spin correlation 
function is exactly the same. The ground state is unique and has antiferromagnetic long- 
range order. 
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